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SUMMARY 

The  simple  line-vortex  model  of  laterally-symmetric  separated  flow 
past  a  circular  cone  at  incidence,  originally  proposed  by  Bryson,  is 
studied  in  detail,  both  numerically  and  -nalytically  through  asymptotic 
expansions.  Two  branches  of  solutions  are  found,  an  'upper'  branch  which 
includes  Bryson's  original  results,  and  a  lower  branch,  on  which  the 
solutions  are  found  to  be  physically  unrealistic.  The  physically  realistic 
branch  is  shown  to  yield  multiple  solutions  for  separation  lines  to 
windward  of  a  critical  position.  A  proposed  modification  to  the  condition 
imposed  on  the  flow  at  the  separation  line  is  shown  to  have  little  effect 
on  the  behaviour  of  the  model. 
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1  INTRODUCTION 

There  is  considerable  interest  in  the  aerodynamics  of  aircraft  and  missiles  at  high 
angle  of  attack.  The  flow  past  the  nose  of  such  a  vehicle  will  often  exhibit  vortex 
separation.  This  separation  usually  starts  at  the  pointed  apex,  and  hence  considerable 
information  might  be  obtained  from  studying  flows  past  cones.  As  these  flows  involve 
separation  one  ought,  in  principle,  to  solve  the  full  Navier-Stokes  equations.  However, 
Fiddes'  has  shown  that  an  inviscid  model  involving  vortex  sheets  can  give  useful  predic¬ 
tions  of  real  flows  involving  separation  from  smooth  surfaces,  when  the  separation  line 
is  determined  by  a  calculation  in  which  the  interaction  with  the  viscous  flow  is 
represented. 

2 

Bryson  introduced  a  very  simple  inviscid  model  for  the  flow  past  a  slender  circular 
cone  at  incidence.  He  made  the  slender-body  assumption  and  modelled  the  separated  flow  by 
isolated  line-vortices.  Each  such  vortex  is  linked  by  a  cut  to  the  separation  line  from 
which,  in  the  real  flow,  it  would  be  fed.  It  is  required  that  the  net  force  on  the  vortex 
and  cut  should  vanish.  The  position  of  the  separation  line  must  be  specified,  and  Bryson 
required  that  this  line  should  be  a  streamline  of  the  three-dimensional  flow. 

Although  this  model  is  too  simple  to  be  of  practical  use,  its  properties  are  of 
interest  because  on  the  one  hand  it  represents  a  degenerate  case  of  the  vortex-sheet  model 
and  on  the  other  hand  its  solutions  may  provide  initial  approximations  for  the  numerical 
procedures  needed  to  calculate  the  behaviour  of  the  vortex-sheet  model.  Its  simplicity 
allows  a  complete  investigation  of  its  solutions. 

The  study  described  in  section  2  reveals  a  'lower'  branch  of  solutions,  in  addition 
to  the  upper  branch  which  contains  Bryson's  original  results.  On  this  lower  branch  the 
vortex  lies  on  the  windward  side  of  the  'separation  line'.  It  is  shown  that,  for  the 
solutions  on  this  branch,  the  postulated  'separation  line'  is  actually  an  attachment  line, 
so  that  these  solutions  must  be  discarded.  It  emerges  that,  even  on  the  physically 
realistic  branch,  multiple  solutions  are  possible  if  the  separation  line  lies  to  windward 
of  a  critical  position,  which  is  itself  to  windward  of  the  position  Bryson  chose  for  his 
numerical  investigation.  Physical  arguments  are  put  forward  for  preferring  one  of  these 
solutions . 

The  condition  that  the  separation  line  should  be  a  streamline  of  the  three- 
...  .  .  .3 

dimensional  flow  is  somewhat  arbitrary.  Smith  has  proposed  an  alternative  condition 
which  is  consistent  with  the  vortex-sheet  model  of  the  flow.  In  section  3,  the  effect  on 
the  solution  of  using  the  alternative  condition  is  shown  to  be  small.  The  main  features 
revealed  in  section  2  are  unchanged. 

2  BRYSON'S  MODEL 

2 

For  completeness,  the  model  of  Bryson  is  re-derived  in  section  2.1,  before  the 
solution  techniques  and  results  are  described  in  the  remainder  of  section  2. 

2. 1  Equations 

Some  of  the  notation  to  be  used  is  illustrated  in  Fig  1.  A  Cartesian  coordinate 
system  is  adopted,  and  the  apex  of  the  cone  is  assumed  to  lie  at  the  origin  of  the 


coordinate  system.  The  axis  of  the  cone  is  taken  to  lie  along  Ox  ,  and  the  coordinates 
in  any  plane  perpendicular  to  this  axis  are  y  and  z  . 

The  slender-body  assumption  is  made.  Hence  the  velocity  in  a  plane  perpendicular 
to  the  axis  of  the  body  (the  cross-flow  plane)  may  be  described  by  a  potential  which 
satisfies  Laplace's  equation  in  the  cross-flow  plane,  and  so  is  the  real  part  of  a  complex 
analytic  function.  The  equations  of  the  model  will  therefore  be  formulated  in  terms  of  a 
complex  variable,  oi  =  y  +  iz  .  In  the  slender-body  description  of  the  attached  flow  past 
a  conical  body,  the  component  of  the  velocity  in  the  cross-flow  plane  is  conical,  in  the 
sense  of  being  constant  along  rays  through  the  apex.  It  will  be  assumed  that  the  separated 
flow  is  also  conical;  in  particular,  that  the  separation  line  is  straight,  and  that  the 
line-vortices  which  represent  the  real  separation  vortices  are  also  rectilinear. 

s  ,  the  radius  of  the  circular  cross-section  of  the  cone  in  the  cross-flow  plane, 

is  equal  to  x  tan  y  ,  where  y  is  the  (small)  semi-angle  of  the  cone.  The  separation 

lines  are  co  =  uj  s  where  u  =  e10  or  -e  10  ,  and  6  ,  which  satisfies  -90°  <  6°  <  90°, 
s  s 

is  termed  the  separation  angle.  wvs  and  -ui^s  (where  a  bar  denotes  the  complex  conju¬ 
gate)  are  the  positions  of  the  vortices  modelling  the  symmetric  separation.  The  velocity 
in  the  far  field  is  (U  cos  a,  0,  U  sin  a)  where  a  is  the  (small)  angle  of  incidence 
in  radians,  and  U  is  the  free  stream  speed.  Because  the  flow  is  conical  the  circulation 
of  the  vortices  is  proportional  to  x  and  is  taken  to  be  equal  to  2irrsU  tan  y  . 

The  condition  that  the  component  of  the  velocity  normal  to  the  surface  of  the  cone 
is  zero  gives 

v^  =  V  tan  y  , 

where  is  the  velocity  component  on  the  circle  along  its  outward  normal  in  the  (y,z) 

plane.  The  complex  potential,  W  ,  can  be  seen  to  be: 


The  first  term  represents  the  effect  of  the  streamwise  growth  of  the  cross-section,  the 

second  term  describes  the  attached  flow  round  the  cross-section  of  the  cone,  and  the 

third  term  gives  the  effects  of  the  two  vortices  and  their  images  in  the  circle.  The 
.  .  .  2  . 

condition  which  Bryson  imposes  on  the  flow  at  the  separation  line  is  that  the  velocity 
there  is  parallel  to  the  separation  line.  In  the  real  viscous  flow,  the  limiting  direc¬ 
tion  of  the  flow,  as  the  wall  is  approached,  is  parallel  to  the  separation  line.  For  an 

3 

inviscid,  vortex-sheet  model  of  the  separated  flow,  Smith  points  out  that  the  flow  on 
the  downstream  side  (but  not  on  the  upstream  side)  of  the  separation  line  must  be  parallel 
to  it.  An  alternative  condition  valid  for  both  vortex-sheet  and  line-vortex  models  is 
explored  in  section  3.  Bryson's  condition  leads  to: 
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A  +  s2\  _  Ts  tan  y  |~  1 _ 1 _ 1  +  1  *| 

H  2 )  i  u-ti)S  w  -  (s/S  )  w  +  sE  u  +  (s  /u  ) 

\ll)/  L  v  v  v  v_| 


at  w  »  wgs  ■  se  .  This  equation  is  satisfied  in  a  trivial  sense  for  to  =  ±is  . 


Other¬ 


wise  it  implies: 


T  tan  y  = 


a(i  +  ^ve~l9)  (>  ('  ~y~ie)  ('  ~yie) 


(co  5  -  1)  (w  +  E  ) 
V  V  V  V 


Bryson  assumes  that  the  line  vortex  is  inclined  at  a  small  angle  to  the  local 
velocity  vector  so  that  the  resulting  force  on  the  vortex  exactly  balances  the  force  on 
the  cut.  For  the  starboard  vortex,  this  leads  to: 

(wvs  -  ugs)  ^  (2nTsU  tan  y)  +  (2nFsU  tan  y)  (uvs)  ”  (2trrs  tan  y) (vv  + iwv)  =  0  ,  (2-4) 

where  the  first  term  comes  from  the  force  on  the  cut,  the  second  term  from  the  force  on 

the  vortex  due  to  the  free-stream  and  the  final  term  from  the  force  on  the  vortex  due  to 

the  cross-flow,  v  + iw  is  the  cross-flow  velocity  at  the  vortex  at  us,  and  is 
v  v  1  v 

given  by 

dW 

v  -  iw  *  — j —  (2-5) 

v  v  du 

U“SU 

V 

where  W,  *  W  -  1  ln(u-u  s)  ,  the  derivative  of  W.  ,  rather  than  W  ,  being 

1  I  V  1 

required  because  a  straight  line-vortex  has  no  self-induced  velocity. 

Substituting  equation  (2-5)  into  the  complex  conjugate  of  equation  (2-4)  and 
recalling  that  s  »  x  tan  y  leads  to 

tan  y(2S  -e_i9)  -  -  iafl  ♦  ±)  *  +  liff.J  l  - >,  -  - -  *  -  '1 

'  V  /  \  U2/  “v  1  L  “v-  (1/V  “v  "v  uv  (1/ulv)J 


so  that 


\  v  w  /  l  2  /  w  ♦  <*>  /  2  \,  -  IX 

N  V/  V  U>  /  V  V  I  b)  +  1 1  (u  (*>  ~  I ) 

V  \  V  /  V  v 


where  a  *  a/tan  y  >  and  is  termed  the  incidence  parameter. 
Substituting  f  from  equation  (2-3)  into  this  gives 
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"(25v-'ie-^)  -  - '('  *4) 


i(t  *5/xe)(i  +  yl9)  (i  -  y~l6)  Q  ~ y'9) 


2 

(u>  0)  -  1)  (u>  +  w  ) 

V  V  V  V 


i(i  +  ^ve~l9) (i  ♦y8l9)(1  ~  v"16)!1  ~y19) 

(«aj*  t)(uvSv“1)2 


.  (2-7) 


0 

If  t  *  u>  e  this  becomes 
v 


a"'e-ie^ 


ie(2i  . ,  _  I)  =  -  i(,  ♦  £^L)  ♦ 

'  /  \  t  /  (tt-l)Cte1  +  te 


-t)(l  -t) 
ie)2 


.  i(1  + te~2i6)(l  +  teZl9)  ( 1  -t)(l  -t)  -2  -2ie. 

(i  +  tVift)(.  +  t2.-2ie)(tt-.)2'' 


(2-8) 


Equations  (2-7)  and  (2-8)  are  each  equivalent  to  two  real  equations  relating  the  incidence 
parameter,  a  ,  and  the  real  and  imaginary  parts  of  the  complex  vortex  position,  or 

t  respectively.  It  is  convenient  to  proceed  by  specifying  the  imaginary  part  of  or 

t  and  then  solving  to  find  the  incidence  parameter  and  the  real  part  of  or  t 

respectively.  From  a  knowledge  of  these  three  parameters,  F  can  be  evaluated. 

The  final  equation  needed  is  an  expression  for  the  normal  force.  Let  be  the 

normal  force  on  the  cone  forward  of  the  plane  x  =  constant,  non-dimensionalised  with 

5. 
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2  2 

the  cross-sectional  area  of  the  cone  (=  tts  )  and  JpU  where  U  is,  as  before,  the  free- 
stream  velocity  and  p  is  the  density  of  the  fluid.  Then  it  can  be  shown  thad 


“N 


tan  y 


=  2a  +  8Ty: 


■(' '  ?) 


(2-9) 


where  m  =  y  +  iz  and  w 
v  v  v  v ' 

tion  are  known. 

2.2  Numerical  methods 


R  .  can  be  found  if  the  vortex  position  and  circula- 


It  is  possible  to  apply  three  criteria  to  assess  the  relevance  of  solutions  of 
equations  (2-7)  and  (2-8).  It  is  assumed  that  the  y  and  z  coordinate  directions  are 
chosen  so  that  a  is  positive.  Hence  one  restriction  must  be 


>  0  . 


(2-1 0a) 


The  vortices  modelling  the  separation  must  lie  outside  the  cone  (and  their  images  inside 
the  cone) ,  and  so 


>  1 


(2- 10b) 


mj. 
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If  the  cut  is  regarded  as  the  vestige  of  a  vortex  sheet,  it  is  necessary  to  ensure  that 
the  two  cuts  joining  the  vortices  to  their  respective  separation  lines  do  not  cross.  For 
symmetric  separation  (in  which  the  two  cuts  are  reflections  of  each  other  in  the  line 
y  =0)  it  is  therefore  necessary  that  the  y  coordinate  of  the  vortex  and  its  associated 
separation  point  have  the  same  sign.  This  leads  to 

<R(wv)  >  0  ,  (2-10c) 

because  -90°  <  0  <  90°  . 

First  the  solution  of  equation  (2-7)  will  be  considered.  It  will  be  assumed  that 

y(u>  )  is  specified.  If  real  and  imaginary  parts  of  equation  (2-7)  are  taken,  two  equa- 

.  V  .  -1  -1 
tions  relating  a  and  (R(uv)  result.  If  a  is  eliminated  between  these  two  equa¬ 
tions  a  !7th-order  polynomial  for  y  =  R(u  )  can  be  obtained.  Taking  z  =  ./(u  )  and 
2  2  2  .  v  .  v  v  v 

R  =  y^  + zv  ,  the  polynomial  can  be  written  as 


2y  z  (- 
v  v\ 


2  2  \  2 
2z  R  +  sm  6  R  +  z  )4y  x 


[■*('  *v'i8)0  V)0  -Vie)  -  (r2-i)2("J*,)('"2v*‘)] 


(2yyR2  -  R2  cos  0  -  y^ 


+  r4(r2~1)(U)v+1)(“v+1)(I  +V'l9)('  +“vel6)('  "V  l9)('  ‘Vi8) 

+  r444('  +vie)('  + ')]  ■ 


(Despite  its  appearance  this  expression  is  real  and  can  be  written  as  a  polynomial  in 
yv  .  It  has  been  left  in  this  form  for  brevity.) 

All  roots  of  this  were  found  using  the  method  of  Grant  and  Hitchins  as  implemented 
•  •  A 

by  the  Numerical  Algorithms  Group  in  their  subroutine  C02AEF.  Unacceptable  solutions 
for  fl(u)v)  were  then  eliminated  using  the  requirement  that  they  must  be  real  (the  poly¬ 
nomial  may  have  some  complex  zeros,  even  although  all  the  coefficients  are  real),  and 
the  criteria  of  equations  (2-10). 

The  method  seems  to  work  well  in  general,  but  was  poor  near  the  surface  of  the  cone. 
This  is  due  to  the  occurrence  of  repeated  roots  of  the  polynomial  at  the  surface,  because 
the  subroutine,  although  second-order  convergent  for  single  roots,  exhibits  numerical 
first-order  convergence  for  repeated  roots. 
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When  the  method  just  described  for  the  solution  of  equation  (2-7)  was  unsatisfactory, 
equation  (2-8)  was  solved  (using  FORTRAN  with  double-precision  real  variables)  for  speci¬ 
fied  J(t)  ,  to  find  just  one  root,  in  the  following  way.  By  taking  real  and  imaginary 
parts  of  equation  (2-8)  two  expressions  for  a  can  be  found.  Subtracting  these  two 
expressions  gives  a  function  ip  ,  say,  of  g  =  fl(t)  ,  the  zeros  of  which  give  the  roots 
of  equation  (2-8).  To  find  a  zero  a  modified  interpolation  scheme  was  used,  in  which  two 
values  of  £  must  be  supplied,  say  £A  and  ,  such  that  the  signs  of  ilKf^)  and 
iKCg)  differ.  Then  there  is  at  least  one  root  of  iJj=0  between  £A  and  and  the 

scheme  finds  one  root.  Since  the  method*  does  not  seem  to  be  well  known  it  is  explained 
in  Fig  2. 

Although  the  scheme  is  similar  to  an  ordinary  interpolation  scheme,  which  does  not 
involve  ordinate  halving  (ie  the  instructions  Fg  =  Fg/2  and  F^  =  would  be  omitted 

from  the  bottom  two  boxes  in  Fig  2b),  experience  has  shown  that  the  modified  scheme  is 
often  much  faster  and  rarely  much  slower  than  the  ordinary  scheme.  It  has  the  additional 
advantage  that  any  root  can  be  found  to  a  known  accuracy,  because  the  iterations  are,  in 
general,  stopped  when  | £A  -  £g (  is  sufficiently  small,  and  the  error  in  the  zero  has  an 
upper  bound  of  |g  -£  |  '  *n  t*ie  OT^naT^  scheme  the  root  is  said  to  have  been  found 
when  tyiO  is  sufficiently  small.  As  |£  -£g|  may  still  be  large  it  is  difficult  to 
estimate  the  accuracy  to  which  the  root  has  been  found. 

2 . 3  Asymptotic  solution  when  the  vortex  is  near  the  surface  of  the  cone 

In  this  section  equation  (2-8)  will  be  studied,  rather  than  equation  (2-7) ,  because 
the  asymptotic  analysis  appears  to  be  slightly  simpler.  A  similar  analysis  could  have 
been  carried  out  on  equation  (2-7) . 

The  small  parameter  is  taken  as  5  =  ./(t)  .  Details  of  the  asymptotic  analysis  are 
given  in  Appendix  A.  Here  only  the  results  are  given. 

<R(t)  =  1  +  A6  +  tat^-S  —  ^  62  +  0(S3)  , 
a->  =  Llinl  +  (2_c°s_e  .  JOA  Bin  e  V  +  0(62)  >  (2-M) 


where  A=l//3  if  6  >  0 
A  =  -1//3  if  6  <  0, 

so  that  A6  >  0  .  Hence 

a  =  — — y-c-  P.jj  +  ~  cot  e)6  +  0(62)]  ,  (2-12) 

T  =  4a  cos  6 (Act)  +  0(62)  ,  (2-13) 


*  The  method  of  solution  was  suggested  to  the  author  by  Dr  F.A.  Johnson,  while  at 
RRE  Malvern. 
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from  equation  (2-3) ,  and 

Q 

- —  =  2a  +  0(52)  ,  (2-14) 

tan  Y 

from  equation  (2-9).  These  results  are  in  good  agreement  with  the  numerical  results  for 
sufficiently  small  } 6 )  . 

For  given  0  ,  the  locus  of  possible  vortex  positions  (a  varying  along  the  locus) 
appears  to  have  two  branches  near  the  surface  of  the  cone.  Their  asymptotic  form, 
to  0(6),  is  shown  by  the  dashed  lines  in  Fig  3,  for  the  two  cases  6  =  56°  (the  case 
studied  by  Bryson2)  and  6  =  14.3°  (=  0.25  radian),  the  actual  loci  being  shown  by  the 
dotted  curves.  Also  noted  on  Fig  3  is  the  behaviour  (as  the  vortex  moves  away  from  the 
surface  of  the  cone)  of  the  incidence  parameter,  the  circulation  and  the  normal  force. 

The  critical  value  of  0,  0^  ,  which  determines  the  change  of  behaviour  of  F  with  a  on 
the  upper  branch,  can  be  found  from  equation  (2-11).  It  is  given  by 


tan  ec  =  (#)(t§a) 


3/3 

5 


as  A  =  +  1//^  on  the  upper  branch,  so  that  0^  =  46.1  .  A  similar  change  of  behaviour 
of  T  with  a  can  be  observed  on  the  lower  branch,  the  critical  0  here  being  -0^  . 

In  order  to  clarify  the  nature  of  the  upper  and  lower  branches  of  the  solution,  it 
is  helpful  to  examine  the  flow  near  the  separation  line.  The  complex  conjugate  of  the 
cross-flow  velocity  is  given  by  the  derivative  of  (2-1): 


.  dW 

v  -  iw  =  -r—  =  U  tan  y 

dm 


[js  .  A  +  s2N  +  Vs  /  1 _ fv _ 1  +  “v  ^1 

lu  \  m2/  ^  (“'“v®  wwv~s  w  +  s5v  umv  +  syl 


(2-15) 
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As  in  (2-8),  the  position  of  the  starboard  vortex  is  written  as  «v  =  te  ,  where,  by 
(2-11) 


t  =  1  +  6^i  +  +  0(62) 


and  the  upper  and  lower  signs  refer  to  the  upper  and  lower  branches,  on  which  6  is  posi¬ 
tive  and  negative,  respectively.  The  velocity  is  required  at  points  which  are  closer  to 
the  separation  line  than  is  the  vortex,  so  we  write 


=  s  ( I  +  ?)e 


10 


where  | C  j  <  j  6  |  . 

With  r  =  (4//3)a|6|  cos  0  ,  from  (2-13),  we  find  the  leading  terms  of  (2-15)  are 


v  -  IV 


U  ta*-  o 


ie(.  ♦  3a 4  ...)  . 


(2-16) 
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At  a  point  in  the  fluid  on  the  normal  to  the  surface  of  the  cone,  t.  is  real  and  positive. 
The  velocity  along  the  outward  normal  to  the  surface  of  the  cone  is  given  by  a  contribu¬ 
tion  (v  cos  8  +  w  sin  8)  cos  y  from  the  cross-flow,  together  with  a  component 
-Usin  v(l  +  0(c))  from  the  free  stream,  ie  by  (2-16): 

3aUC  sin  y  cos  6 
6 


to  leading  order  in  c/6  .  The  flow  is  therefore  away  from  the  surface  for  6  >  0  and 
towards  the  surface  for  6  <0  .  Hence  the  'separation  line'  can  represent  a  genuine 
separation  line  for  the  upper  branch  of  the  solution,  but  represents  an  attachment  line 
for  the  lower  branch  of  the  solution.  The  lower  branch  must  therefore  be  discarded  as 
physically  irrelevant. 

2.4  Numerical  results 

It  has  already  been  noted  that  two  branches  of  the  locus  of  vortex  positions  for 
fixed  0  ,  only  the  upper  one  of  which  corresponds  to  separation  at  0  ,  are  obtained 
from  the  asymptotic  solution  when  the  vortex  is  close  to  the  cone.  The  two  branches 
obtained  from  the  solution  of  equation  (2-7)  by  the  method  of  section  2.2  are  shown  by  the 
dashed  curves  in  Fig  4  for  0  =56°  (the  case  studied  by  Bryson^).  Whilst  only  the  upper 
branch  is  physically  relevant,  for  completeness  results  for  both  branches  are  shown  in 
Fig  4  and  in  subsequent  figures. 


By  solving  either  equation  (2-7)  or  equation  (2-8),  Figs  5  and  6  can  be  obtained. 
These  two  figures  show  loci  of  possible  vortex  positions  with  either  0  varying  and 
a  constant,  or  a  varying  and  0  constant.  In  Fig  5  attention  is  restricted  to  the 
upper  branches  (at  constant  8),  while  in  Fig  6  attention  is  restricted  to  the  lower 
branches.  All  the  solutions  for  a  >  0  lie  in  the  upper  quadrant,  as  shown.  If  a  <  0  , 


the  mirror  image  in  the  horizontal  axis  of  Fig  5  or  6  is  obtained.  If  a  locus  for 
constant  0  ,  say  0  =  -14.3°  ,  in  Fig  5  is  followed  into  the  lower  quadrant,  it  con¬ 
tinues  smoothly  as  the  mirror  image  of  the  locus  for  9  =  14.3°  in  Fig  6,  ie  as  a  locus 
of  lower  branch  solutions.  Conversely,  the  lower  branch  loci  of  Fig  6  continue  as  the 
mirror  images  of  the  upper  branch  loci  of  Fig  5. 

The  possible  non-uniqueness  of  solutions  of  the  Bryson  model  is  clearly  illustrated. 
If  attention  is  restricted  to  the  lower  branches,  Fig  6,  then,  for  given  0  >0  ,  there 
is  a  unique  solution  for  the  vortex  position  if  a  cosec  9  ,  and  no  solution  if 
a  cosec  0  .  However,  for  the  same  9  and  a  (> $  cosec  9),  there  will  be  at  least 
one  upper  branch  solution.  Fig  5,  for  the  vortex  position,  and  all  the  upper  branch  solu¬ 
tions  will  differ  from  the  lower  branch  solution.  For  the  upper  branch  (Fig  5),  corre¬ 
sponding  to  physically  realistic  solutions,  the  existence  of  solutions  now  depends  on 
9  and  a  in  a  more  complicated  way,  which  is  illustrated  in  Fig  7.  A  curve  a  =  a^(6) 
represents  the  lowest  value  of  a  for  which  a  solution  with  separation  at  9  exists. 

If  9  >  9c  —  46.1°  ,  then  a^(9)  -  $  cosec  9  ,  the  same  boundary  as  was  found  for  the 
existence  of  lower  branch  solutions,  and  there  is  \  unique  solution  for  a  >  a^(9)  » 
with  no  solution  for  a  <  a^(9)  •  If  0  <  9  <  0^,  a^(0)  is  a  numerically-defined 
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function,  with  a^(e)  <4  cosec  8  ,  and  there  are  three  possibilities:  for  a  <  a^(6)  , 
there  is  no  solution;  for  a^O)  <  a  cosec  8  ,  there  are  two  solutions,  which  coincide 
for  a  =  a^(8) ;  and  for  a  cosec  8  ,  there  is  just  one  solution.  If  8  <  0  ,  there  is 
no  solution  for  a  <  aL(8)  and  there  are  two  solutions  for  a  >  a^(8)  ,  coinciding  when 
a  =  a^(8)  . 

Figs  5  and  6  also  illustrate  the  regions  in  which  the  vortex  can  lie.  In  Fig  5  the 
locus  of  vortex  positions  for  8  =  -90°  gives  one  of  the  boundaries,  because  of  the 
requirement  that  the  two  (symmetric)  cuts  must  not  cross.  This  locus,  and  the  upper 
branches  of  the  loci  of  vortex  positions  for  fixed  8  with  -90°  <  8  <  0°  are  all 
asymptotically  inclined  at  ir/3  to  the  horizontal.  Their  approach  to  this  asymptote  is 
described  in  Appendix  B.  Another  boundary  is  formed  by  the  Foppl  curve.  This  is  the 
locus  of  points  at  which  vortices  can  lie  symmetrically  in  the  wake  of  a  circular  cylinder 
in  plane  flow  (see  page  223  of  Ref  5,  for  instance),  and  so  corresponds  naturally  to  the 
limit  of  a  -*•  °°  .  Its  equation  is: 


/l  +  y2 


=  1  +  y  +  2y 


so  it,  too,  is  asymptotically  inclined  at  ir/3  to  the  horizontal.  The  remaining 

boundary  (apart  from  the  surface  of  the  cone)  consists  of  part  of  the  line  z  =  0  .  The 

intersection  of  this  line  with  the  locus  of  vortex  positions  for  8  =  -90°  may  be  found 

i  0 

by  substituting  u  =  y  and  e  =  -i  into  equation  (2-7).  Taking  the  real  part  of  the 

.  .  V  .  -1  -1 
resulting  equation  gives  a  =  0  .  Taking  the  imaginary  part,  with  a  =  0  ,  gives  a 

.  .  .  .  2 
quartic  in  y  ,  which  is  a  quadratic  in  y  .  This  may  be  solved  to  give 


•  (« *  ©') 


=  1.878 


In  Fig  6  the  lower  branch  of  the  locus  of  the  vortex  positions  when  8  =  90°  and 
the  line  z  =  0  give  two  of  the  boundaries.  By  the  same  method  as  above  these  can  be 
shown  to  intersect  at 

y  =  1 .878  z  =  0  . 

In  Figs  8  and  9  the  variation  of  the  circulation  parameter  r  with  the  incidence 
parameter,  a  ,  at  constant  separation  position,  8  ,  is  illustrated.  The  abscissa 
starts  at  a  =  I .5  because  there  is  no  solution  to  the  Bryson  model  when  a  <  1.5  . 

In  Fig  8  the  upper  branches  are  considered,  while  in  Fig  9  attention  is  restricted  to  the 

lower  branches.  The  change  in  character  of  the  curves  in  Fig  8  occurs  at  8  =  8c  —  46,1°. 

For  0  <  8  <  8^  ,  there  are  two  values  of  T  ,  corresponding  to  the  two  solutions 

indicated  in  Fig  7,  for  a^(  8)  <  a  <  $  cosec  8  ,  eg  for  4.17  <  a  <6.06  when  8  =  14.3°. 

For  8  >  ec,  F  is  uniquely  determined  (on  the  upper  branch)  by  8  and  a  . 

3  MODEL  WITH  MODIFIED  BOUNDARY  CONDITION  ON  THE  SEPARATION  LINE 

3. I  Equations 

Bryson's  model  requires  the  separation  line  to  be  a  streamline  of  the  three- 
dimensional  flow,  so  that  vt^  ,  the  (mean)  tangential  cross-flow  velocity  component  at 

k . , .  .  . _J 


the  separation  line,  is  zero.  However  by  considering  the  behaviour  of  a  vortex-sheet 

3 

model  at  the  separation  line.  Smith  has  derived  a  condition  on  v^  ,  which  can  be  applied 
whether  the  sheet  is  present  or  not.  For  a  body  of  revolution,  such  as  a  circular  cone, 
with  the  separation  line  along  a  meridian,  Smith's  condition  becomes 


(2irrsU  tan  y) 


in  the  direction  of  increasing  6  ,  using  the  same  notation  as  in  section  2.1. 

If  this  is  the  only  modification  made  to  the  model,  it  may  be  seen  from  the 
potential  given  in  equation  (2-1)  that,  at  the  separation  point  in  the  y-z  plane, 

ie 

uj  »  u  s  =  se  , 

s  * 

(\  .  s^\  „  .  FsU  tan  "V  I  1  1  1  1  ”|  „  r-s-  •  -i( 

"tall  +  — ~  I  U  +  - - - L - 1 — - - — sr—  +  - 7 — ; — r-  =  ~  U  tan  y/irT  ie 

\  2/  l  (jj-uis  u  -  (s/tii  )  u  +  sw  u  +  (s/w  )  I 

'01/  L  v  v  v  v_J 

instead  of  equation  (2-2),  where  the  positive  square  root  is  to  be  taken.  This  gives, 
instead  of  equation  (2-3), 


r  (u  o)  - 1 )  2  cos  e(oiv  +  Sv) 


_  -iew.  .  Te\7  -ie\7.  I  TeV  ~  2a  cos  0  +  ^  =  0 
wve  JV +V  ‘V  )V~\e  ) 


Equation  (2-6),  obtained  from  the  requirement  that  the  net  force  on  the  vortex  and  cut  is 
zero,  is  unaltered. 

—  £  Q 

If  t  =  ,  equation  (3-1)  becomes 


T(tt-l)  2  cos  8(te*~9+te  ^9) 

(1  +  te“2l9)(l  +  te2l9)(l  -t)(l  -t) 


-  2a  cos  6  +  /irF  =  0 


and  equation  (2-6)  may  be  written  in  the  form 


-  J,  +  e  2l6'\  +  ir _  ir(tei9  +  te~l9) 

\  t 2  /  (te18 +  te~iB)  (tt- l)(l  +t2e2i9 


The  normal  force  is  still  given  by  equation  (2-9). 

If  the  imaginary  part  of  t  is  specified,  equations  (3-2)  and  (3-3)  may  be  solved 
to  find  E,  a  and  the  real  part  of  t  .  From  a  knowledge  of  these  the  normal  force  may 
be  calculated. 

3.2  Numerical  methods 

The  method  of  solution  of  equations  (3-2)  and  (3-3)  is  similar  to  that  used  to 
solve  equation  (2-8). 
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The  three  criteria  of  relevance  (equations  (2-10)  of  section  2.2)  still  apply  and 
for  the  same  reasons.  It  is  convenient  to  specify  J(t)  . 

For  any  guessed  <l(t),  r/a  is  evaluated  from  the  imaginary  part  of  equation  (3-3) 
and  then  a  determined  from  the  real  part  of  this  equation.  These  values  are  substituted 
into  the  expression  on  the  left-hand  side  of  equation  (3-2),  taking  the  positive  square 
root.  The  modified  interpolation  scheme  described  in  section  2.2  and  Fig  2  is  then 
employed  to  find  <?(t)  so  that  equation  (3-2)  is  satisfied. 

This  method  does  not  have  the  advantage  of  ensuring  that  all  possible  solutions  of 
the  equations  are  found,  but  a  thorough  search  failed  to  reveal  any  relevant  solutions 
other  than  those  referred  to  in  sections  3.3  and  3.4. 


As  in  section  2.3,  the  small  parameter  is  6  =  ./(t)  .  Details  of  the  asymptotic 
analysis  are  given  in  Appendix  C.  Here  the  results  are  given. 


«(t)  =  1  +  A6  - 


2ttA6  tan  8 


'  8A  tan  8 


4n  tan  8 
81 


where  A  =  +1//3  if  6  >  0, 
A  =  -1//3  if  6  <  0, 

so  that  A6  >  0  .  Hence 


and 


a 


tan  y 


[i  ♦  -  cot  e)«  -  o(fi*)] 

2a  +  0(«2)  . 


r-^ 

O' 


o 

u 

a> 

< 

£ 


These  results  (which  are  only  applicable  when  8  >  0  because  of  the  restriction 
a  >  0)  agree  to  0(6)  with  the  results  given  in  section  2.3,  but  the  next  term  is  of 
different  form. 


For  given  8  the  locus  of  possible  vortex  positions  (a  varying  along  the  locus) 
again  has  two  branches  near  the  surface  of  the  cone. 


The  leading  term  of  the  solution  is  again  shown  by  the  dashed  lines  in  Fig  3,  but 
the  dotted  curves  no  longer  give  the  actual  loci  of  the  vortex.  The  critical  value  of 


is  still 
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The  previous  analysis  of  the  velocity  near  the  point  m  =  se  holds,  and  it  follows  that 
the  lower  branch  solutions  must  again  be  discarded. 

3.4  Numerical  results 

Again  it  was  found  that,  because  of  the  restriction  a  >  0  ,  all  solutions  satisfied 
./(u>v)  >  0  .  For  constant  0  the  locus  of  possible  vortex  positions  still  has  two 
branches,  as  is  illustrated  by  the  solid  curves  in  Fig  4,  for  the  case  6  =  56°  .  Fig  4 
also  illustrates  the  effect  of  modifying  the  condition  at  the  separation  line. 

In  Fig  10  the  variation  of  the  circulation  parameter,  T  ,  with  the  incidence  para¬ 
meter,  a  ,  for  0  =  56°  is  illustrated  by  the  solid  curves.  This  figure  also  shows  the 
effect  of  modifying  the  condition  at  the  separation  line. 

As  the  effect  of  modifying  the  boundary  condition  is  quantitative  rather  than  quali¬ 
tative  all  the  comments  of  section  2.4  are  still  valid.  The  critical  0,  0C  ,  is  still 
tan  '(3/3/5)  and  it  is  still  possible  to  have  non-unique  solutions. 

6  3 

In  his  example  III.l,  Nangia  also  considers  Smith's  boundary  condition.  Fig  4 
here  and  Fig  8  of  Nangia  agree  to  within  the  accuracy  of  curve  plotting. 

4  DISCUSSION 

Bryson's  model  for  separated  flow  round  a  circular  cone  has  been  investigated,  and 
the  occurrence  of  multiple  solutions  clarified.  For  fixed  separation  position,  0  , 
a  branch  of  the  locus  of  the  vortex  position  not  previously  reported  has  been  found.  How¬ 
ever  it  has  been  shown  that  for  this  branch  (termed  the  'lower'  branch  in  this  paper)  the 
so-called  separation  line  is  in  fact  an  attachment  line.  Hence  this  branch  of  the  solu¬ 
tion  can  be  discarded  on  physical  grounds. 

Then,  for  0  >  0^  (0^  =  tan  '(3/3/5)  =  46.1°),  there  is  a  unique  solution  for  any 

given  0  and  incidence  parameter,  a  ,  provided  a  cosec  0  .  If  a  <$  cosec  6 
there  is  no  solution.  The  incidence  parameter  and  the  strength  of  the  vortex  increase  as 
the  position  of  the  vortex  moves  away  from  the  surface  of  the  cone.  It  is  for  a  value  of 
0  (56°)  in  this  range  that  Bryson^  calculated  solutions  of  his  model,  finding  the 
realistic,  upper  branch  of  the  present  investigation.  For  0  <  0  <  0^  with 
a  (0)  <  a  <  i  cosec  0  ,  and  for  -90  <  0  <  0°  with  a  (0)  <  a  there  are  two  possible 
vortex  positions  and  strengths,  so,  as  a  increases  from  a^(0)  with  0  fixed,  there 
are  two  branches  of  the  solution.  Along  one  of  them  the  vortex  moves  away  from  the  cone 
and  becomes  stronger  as  a  increases.  If  0  <  0  <  @c  ,  this  branch  is  a  smooth  continua¬ 
tion  of  the  unique  solution  which  exists  for  a  >  i  cosec  8  .  It  approaches  the  Foppl 
curve  as  a  -*  ■»  .  Along  the  other  branch  the  vortex  moves  towards  the  cone  and  becomes 
weaker  as  a  increases.  This  is  physically  unrealistic.  For  0  >0  ,  the  surface  of 
the  cone  is  reached  when  a  =  f  cosec  0  . 

In  a  real  flow,  the  separation  position,  0  ,  must  depend  on  the  incidence  para¬ 
meter,  a  ,  through  the  viscous  mechanism  governing  separation.  It  is  likely  that,  as 
a  increases,  separation  first  occurs  near  0  »  n/2  ,  and  that  0  falls  as  a  increases. 
By  the  time  8  has  fallen  as  far  as  the  vortex  will  be  well  away  from  the  surface, 

so  that,  when  8  falls  below  0c  ,  the  vortex  will  be  on  the  physically  realistic  branch 
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referred  to  above.  In  terms  of  the  bifurcation  locus  a^(8)  and  t*le  separation  position 
8(a)  in  the  viscous  flow,  it  is  expected  that  a^(8(a))  is  appreciably  less  than  a 
itself,  so  that  there  is  no  way  in  which  the  non-physical  branch  can  be  reached. 

3 

Smith's  modified  condition  on  the  separation  line  has  been  shown  to  have  a  small 
quantitative,  but  no  qualitative,  effect  on  the  results.  Hence  this  modification  has  no 
effect  on  any  of  the  comments  made  in  this  section. 
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Appendix  A 

ASYMPTOTIC  SOLUTION  OF  BRYSON  MODEL  WHEN  THE  VORTEX 
IS  NEAR  THE  SURFACE  OF  THE  CONE 


The  equation  to  be  solved  is 


-1 


/  -2i8\ 

(it  -  1  -  -) 

s  -  \ 

l(i  +  - — —  ) 

\  t/ 

V  t 2  / 

r  -2i8, 


210, 


;i  -  o 


id  +  te~2l6)(l  +te2l6)(l  -t)(l  -t)(l  uV2l9l 
(1  +  t2e2l6)(l  +  t2e~2l0)(tt  -  l)2 


(A-l) 


which  is  equation  (2-8),  given  in  section  2.1. 

Let  the  small  parameter  be  6  =  and  let  fl(t)  =  1  + q  where  q  is  small.  It 

can  be  shown  that  acceptable  solutions  can  only  be  obtained  when  q  =  0(6)  .  Hence  let 

q  -  A6  +  0(62) 

where  A  is  a  constant  to  be  determined.  To  lowest  order  the  left-hand  side  of  equation 
(A-l)  is  then 

a  'e  ( 3A  —  i )  6 

and  the  right-hand  side  is 

-  i()  +e~2i6)  4-  y^2.Vl  u+e-2ie)  . 

4A^ 

2 

By  comparing  these,  it  may  be  seen  that,  since  the  left-hand  side  is  0(6),  A  ■  1/3  .  As 
| t | 2  =  I  +  2A6  +  0(62)  and  |t|  >  1  , 


where  A  *  +1//3  if  6  >  0, 
and  A  =»  -I//3  if  6  <0. 

Now  let 


fl(t) 


1  +  A6  +  B62  +  0(63)  , 


i  0 

multiply  equation  (A-l)  by  e  ,  and  expand  it  to  order  6  .  Taking  real  parts  gives 


-1  2  sin  0  .  ,, 

a  -  - z -  +  0(6) 


Then,  by  considering  the  imaginary  part,  it  may  be  shown  that 

„  _  (8A  tan  6-4) 

B  - £ -  . 


tticZLui^a  I,.,*  bjukij-wr  noma 
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Appendix  A 


Now,  writing 

2  »f-  6  +  DS  *  0(62) 

2 

times  equation  (A-l)  to  order  S  it  is  found  that 

2  cos  6  _  IQA  sin  6 
3  ~  9 

Thus  the  asymptotic  expansions  given  in  section  2,3  are  obtained,  equations  (2-3)  and  (2-9) 
of  section  2.2  being  used  to  obtain  the  asymptotic  forms  of  r  and  . 


-1 

a 

10 

and  expanding  the  real  part  of  e 
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Appendix  B 


ASYMPTOTIC  FORM  OF  THE  UPPER  BRANCH  OF  THE  LOCUS  OF  VORTEX  POSITIONS 
WHEN  THE  SEPARATION  ANGLE,  (TI  IS  ,  FAR  FROM  THE  CONE 


Equation  (2-7),  given  in  section  2.1,  is 


,  \  i(>  +“ve"l6)(1  +V16)('  ~V~19) ('  ~V16) 


(to  n>  —  I )  (to  +  ijj  )  ‘ 
V  V  V  V 


i(1  +  v"19)  0  +yie)('  'V~l6)('  ~v19) 

+  l\  (ui  U)  -  1  )  ^ 

\  V  /  V  V 


If  |«  I  is  assumed  to  be  large  the  left-hand  side  of  this  equation,  to  lowest  order,  is 


and  the  right-hand  side  is,  to  lowest  order. 


\(uv  +  “v)  / 


As  the  right-hand  side  is,  to  this  order,  purely  imaginary 


a  u>  *  o(  1 ) 
v 


and  w  *  pe  to  lowest  order,  where  p  ,  the  modulus  of  w  ,  cannot  be  determined 

v  v 

at  this  order.  Let  1/p  be  the  small  parameter  of  the  expansion.  If 

l>/3  .  Be^.  .  n,  -2."| 

"v  »[*  *  —  *0(0  )J 

and  equation  (B-l)  is  expanded  to  order  o'  it  is  found  that,  to  this  order,  the  right- 
hand  side  of  equation  (B-l)  is  purely  imaginary.  Thus 

-1-  .  -I, 

a  ojv  =  o(p  )  , 

and  p  cannot  be  determined  at  this  order.  ij>  is  an  arbitrary  angle  and  will  be  taken 
equal  to  5 tt/6  as  the  vectors  represented  by  and  e^11^  are  orthogonal.  This 

leads  to 

B  ”  sin  6 


Appendix  B 


Now  let 


u  *  p 

V 


£eW3(l  ♦  +  4  e1*  +  0(p"3)j 


-2  .  .  . 

and  expand  equation  (B— 1)  to  order  p  .As  the  right-hand  side  is  purely  imaginary 


a  'iHv  ”  o(p  3) 


Thus  p  cannot  be  determined  at  this  order  and  ij/  is  an  arbitrary  angle.  It  will,  as 
above,  be  taken  equal  to  5ir/6  .  This  gives 


75  * 


0)  =  p 

V 


[eilt/3(.  +  -  -±3) 


/3pV  p- 


Celx  -4 

♦  +  0(p  ) 


and  expand  equation  (B— 1)  to  order  p  .  If  9=0  it  is  found  that  C  =  0  ,  but  if 
6  <  0  the  real  part  of  the  resulting  equation  gives 


and  hence 


‘  o  •  -  - 

a  =  -  3  sm  9  p  , 


p  =  (-3a  sin  9)^  . 


The  imaginary  part  of  the  equation  gives  a  relationship  between  C  and  x  »  '3U*-  it  is 
not  possible  to  determine  C  and  x  completely  at  this  order. 


For  the  case  9=0,  let 


wv  =  p 


and  expand  equation  (B-l)  to  order  p~  .  Then  the  real  part  of  the  resulting  equation 
gives 

-1  2/3 


and  hence 


p  =  (2^)*  . 


The  imaginary  part  of  the  equation  will  give  a  relationship  between  D  and  y  ,  but  it 
will  not  be  possible  to  determine  D  and  y  completely  at  this  order. 
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Appendix  C 


ASYMPTOTIC  SOLUTION  OF  THE  MODEL  WITH  MODIFIED  BOUNDARY  CONDITIONS 
WHEN  THE  VORTEX  IS  NEAR  THE  SURFACE  OF  THE  CONE 


The  equaCions  to  be  solved  are 


X(tt-l)  2  cos  Site’'®  +  te  *0)  „  „  .  /  ,  -1 

— - h-2 - ~  2  cos  6  +  /ffXa  =  0 

(1  +  te  /10)(1  +  te^1;  ( 1  -  t)  (1  -  t) 


+  _ Xi _  +  Xi(te^6  +  te  *0) 

te10  + te~l6  (1  + t2e2lQ) (tt  -  1) 


which  are  equations  (3-2)  and  (3-3),  given  in  section  3.1,  with  r  replaced  by  Xa  for 
convenience . 

Let  the  small  parameter  be  6  =  t)  and  let  (R(t)  =  1  +q  where  q  is  small. 

It  can  be  shown  that  acceptable  solutions  can  only  be  obtained  when  q  and  X  are  0(d) 
Let 

q  =  Ad  +  higher  order  terms 

where  A  is  a  constant  to  be  determined.  To  lowest  order,  equation  (C-l)  gives 

X  =  (A2  +  1 )  cos  9  4"  • 

To  lowest  order  the  left-hand  side  of  equation  (C-2)  is 

a  1  e  ^0  (3A  -  i)  d  , 


and  the  right-hand  side  is 


.  -2i6^  ,  i(A2  +1)  „  -ie 

-  i(l  +  e  )  +  - = — —  cos  9  e 

2A 


By  comparing  these,  it  may  be  seen  that 


and  hence 


a  (3A  -  i) d  =  0(d) 


.2  .  1 


As  in  Appendix  A,  the  requirement  that  jt|  >  1  gives 


A  *  + 


if  d  >  0 


if  d  <  0 
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X  ■  4A  cos  6  6 


to  lowest  order. 


Inspection  of  equation  (C-l)  now  suggests 


X  -  4A  cos  0  6  ♦  C62  +  0(6  ) 


♦  C6^  +  0(62) 


provided  that  a  1  is  0(1)  .  This  suggests 


4  2 

1  +  A6  +  B67  +  0(6  ) 


Equations  (C-l)  and  (C-2)  expanded  to  0(6*)  give 


A  m  _  A  4nAG6 

3\  cos  6  * 


C  «  4B  cos  6 


where  a  1  »  G  +  0(6^)  and  G  cannot  be  determined  at  this  order.  Letting 

a"'  ■  G  H6^  +  0(6) 


and  expanding  the  real  part  of  equation  (C-2)  to  order  67  gives 


2  sin  6 


H  -  0 


so  that  the  assumption  that  a  is  0(1)  is  justified. 


Now  let 


1  +  A6  «-  B6^  +  E62  +  0(6^) 

4A  cos  0  6  +  C6^  +  F62  +  o(fi^) 


~l  _  2  sin  6 


+  K6  +  0(6^)  . 


From  equation  (C-l)  and  the  imaginary  part  of  equation  (C-2)  to  order  6 


8a  tan  0  4n  tan  0  4 

5  +  81  "  9 


16tt  sin  S  1 6 A  sin  0  16  cos  0 


Then  from  the  real  part  of  equation  (C-2)  to  order  6^ 


2  cos  0  10A  sin  6 
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Thus  the  asymptotic  expansions  given  in  section  3.3  are  obtained,  where  the  equation 


has  been  used  to  find  r  ,  and  equation  (2-9)  has  been  used  for  the  normal  force. 
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LIST  OF  SYMBOLS 


(sgn  6)/3J 


-  a/ tan  y  ,  incidence  parameter 

lowest  value  of  a  for  which  a  solution  exists  for  given  8 
coefficient  of  normal  force,  based  on  cross-sectional  area 


radius  of  circular  cross-section  of  cone 


v  ,w 
v  v 


y  , 2 

•’v  V 


free-stream  speed 

components  of  cross-flow  velocity 

component  of  cross-flow  velocity  normal  to  section  of  cone 
component  of  cross-flow  velocity  tangential  to  section  of  cone 
components  of  cross-flow  velocity  at  position  of  starboard  vortex 
complex  potential 
modified  complex  potential 
Cartesian  coordinates,  see  Fig  ’ 

coordinates  of  starboard  vortex  in  cross-flow  plane 
angle  of  incidence 

(circulation  of  starboard  vortex) /2irUs  tan  y 
semi-angle  of  cone 

=  Jit)  ,  small  parameter  in  expansions 
small  complex  quantity,  section  2.3 
position  of  separation  line,  see  Fig  l 
critical  value  of  8  (tan  *3^ /5) 
argument  of  d* 

-  r/a 

density  of  fluid,  except  in  Appendix  C 
function  of  which  a  zero  is  sought 
•  y  *  iz  ,  complex  variable  in  cross-flow  plane 
value  of  i^/ s  on  starboard  separation  line 


value  of  u/s  at  starboard  vortex 
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Fig  3  Vortex  positions  near  the  cone,  Bryson  model 
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Vortex  positions  for  e=56°,  both  model 


Fig  8  Variation  of  circulation  parameter  with  i 
Bryson  model,  upper  branch 


Pig  9  Variation  of  circulation  parameter  with  incidence  parameter, 
Bryson  model,  lower  branch 
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